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Abstract. We develop a self-consistent model of turbulence in a local inter- 
stellar medium (ISM) . The model describes a partially ionized magnetofiuid ISM 
in which a neutral hydrogen fluid interacts with a plasma dominantly through 
a charge exchange. The ISM turbulent correlation scales in our model are much 
bigger than the shock characteristic length-scales. Unlike small length-scale lin- 
ear collisional dissipation in the fluid, the charge exchange processes can be 
QJ ' effective unpredictably on a variety of ISM length-scales depending upon the 

neutral and plasma densities, the charge exchange cross section and the charac- 
2- teristic length scales. We find, from scaling arguments that the charge exchange 

^2 ' interactions modify spectral transfer associated with large-scale energy contain- 

er , ing eddies. Consequently, the ISM turbulent cascade are steeper than those 

^j • predicted by Kolmogorov's phenomenology. 
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'"5 ■ 1. Introduction and Model Equations 

Small scale turbulence in the local interstellar medium (ISM) is a largely unex- 
plored field given the complexity associated with the ISM turbulent processes 
( Shaikh et. al. 1120061 ). The local ISM comprises of partially ionized, magnetized 
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plasma protons and almost equal number of neutral particles. The plasma and 

the neutral particles in ISM interact mutually through charge exchange. The 
physics of these small-scale turbulent motions is far more complex than ever 
thought. Not only that it holds the key to our crucial understanding of the 
global heliospheric i nteractions such as nature and c haracteristic of heliospheric 
shocks, heating etc ( Zank I [l999l ; iPauls et. aMll995l ). it is also increasingly be- 
lieved to be pivotal to many puzzles of astrophysics including origin and trans- 
port of cosmic rays, Fermi acceleration, gamma-ray bursts, ISM density spectra 
etc. Yet, there exists no self-consistent simulation model that unravels multi- 
component and multiple-scale ISM turbulent phenomena. A prime goal of this 
paper is therefore to develop a self-consistent plasma-neutral ISM turbulence 
model based on analytic methods and numerical simulations. 

The underlying model is based on the following assumptions. Fluctuations 
in the plasma and the neutral fluids are sufficiently isotropic, homogeneous, 
thermally equilibrated and turbulent. No mean magnetic field and velocity flows 
are present at the outset. There however may generate local mean flows due 
to self-consistently excited nonlinear instabilities. The characteristic turbulent 
correlation length-scales are typically smaller than charge-exchange mean free 
path lengths in the ISM flows. Nevertheless, they are large enough to treat 
any localized shocks as discontinuities. In other words, the characteristic shock 
length-scales are relatively small compared to the ISM turbulent fluctuation 
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length-scales, and finally boundary conditions are periodic. It should be further 
noted that the neutrals coming from the solar wind are not considered because 
they tend to anisotropize the distribution functions substantially. Our model 
thus simulates a localized ISM. The fluid model describing nonlinear turbulent 
processes in the interstellar medium, in the presence of charge exchanges forces, 
can then be cast into plasma density (p p ), velocity (U p ), magnetic field (B), 
pressure (P p ), as follows. 



S- + V • ( Pp \J p ) = 0, (1) 



p p (JL + U p • v) U p = -VPp + Ij x B + Q M (U P , V n ), (2) 

/9B 

_ = Vx(UpXB), (3) 

de ( 1 t P c \ 

^ + V- f- / 9 P ^p 2 Up + ^ T ^ /0 pUp + — ExBJ =Q E (V p ,V n ). (4) 

where e = l/2p p Up + P p /(^f — 1) + B 2 /8ir. The above set of plasma equations 
is coupled self-consistently to the ISM neutral density (p n ), velocity (V n ) and 
pressure (P n ) through a set of hydrodynamic fluid equations as below. 

(5) 

r n,Up), (6) 





^ + V • (p n V n ) = 0, 




Pn (g- t + Y n- v) V n = - V P n + Q^ 


d 

at 1 


(I T/2 , Pn \ lV7 (I T/2w 7 Pn 

„PnV n + I + V • p n V n V n + 

\2 7-1/ \2 J-lp n 



PnVn) =Q E (Vn,V p ). (7) 

Equations ([1]) to ([7]) form an entirely self-consistent description of the coupled 
ISM plasma-neutral turbulent fluid. The charge-exchange momentum sources in 
the plasma and the neutral fluids are described respectively by terms Qm(U p , V n ) 
and QA/(V„,Up). Similarly, charge exchange energy sources in the plasma and 
the neutral fluids are given by Qs(Up,V n ) and Q_B(V n ,Up) respectively. In 
the absence of charge exchange interactions, the plasma and the neutral fluid 
are de-coupled trivially and behave as ideal fluids. While the charge-exchange 
interactions modify instantaneous momentum and the energy of plasma and the 
neutral fluids, they tend to conserve density in both the fluids. Nonetheless, the 
volume integrated energy (/ Qe(U p , V n )dv = — J Q_E(V" n , \J p )dv) and the den- 
sity (/ p p dv = const, J p n dv = const) of the entire coupled system will remain 
conserved in the absence of external sources/sinks. 

The ISM turbulence model is normalized using the typical ISM scale-length 
(£o), density (po) and velocity (vq). The normalized plasma density, velocity, 
energy and the magnetic field are respectively; p p = p p /po,\J p = XJ p /vo,P p = 
Pp/poVo,& = B/^oy^TJo. The corresponding neutral fluid quantities are p n = 
p n /pO)U n = XJ n /vo,P n = Pn/poVQ- The momentum and the energy charge- 
exchange forces, in the normalized form, are respectively Q m = Qm^o/Vo^O) Qz = 
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Qe^o/poVQ- The non-dimensional temporal and spatial length-scales are t = 
tvo/£o,~x. = x/£q. The charge-exchange cross-section parameter (a), not ap- 
peared directly in the above set of equations, is normalized as a = n^t^a. 
We define an intrinsic mode in ISM turbulence, the charge exchange mode 
k C e ~ (^o^o )• The latter is different from the characteristic turbulent mode k, 
and is excited naturally when plasma and neu t ral flu ids are coupled in the ISM 
by charge exchange forces [see lShaikh et. aL\ ( 20061 ) for more detail]. 



2. Nonlinear Simulation 

Two-dimensional (2D) nonlinear spectral fluid code is developed to numerically 
integrate Eqs. dU to ([?]). The 2D simulations are not only computationally 
simpler and less expensive (compared with the full 3D), they offer significantly 
higher resolutions even on moderate size small cluster machines like Beowulf. 
The spatial discretization in our code uses discrete Fourier representation of 
turbulent fluctuations based on a spectral method, /(x, t) = X)k/(^^) e_4kX 5 
where k = k x e x + k y e y is a two-dimensional wave vector. The nonlinear decon- 
volution of Fourier modes is performed by computing nonlinear triad interaction 
/(x, t)<7(x, t) = Z)k=k'-k" /(k', *)<7(k", t)S(\i.' — k") which survive for only those 
mode coupling interactions which satisfy Fourier diad constraint k = k' — k". 
These nonlinear interactions in Fourier space conserve rugged invariants of sys- 
tem of Eqs. (pQ) to ([7]). The temporal integration is performed by Runge Kutta 
4 method. The fluctuations in plasma and neutral are initialized isotropically 
(no mean fields are assumed) with random phases and with identical amplitudes 
in the Fourier space. This algorithm ensures conservation of total energy and 
mean fluid density per unit time in the absence of charge exchange and external 
random forcing. Additionally, it satifies the condition of incompressibility asso- 
ciated with the magnetic field, i.e. / |V • B| 2 <iv ~ 10~ 15 at each time step. We 
make use of an ar t ificial scalar potential to achieve V • B = as described in 
Kulikovskii et aL\ (J200ll ). Numerical resolution is 512 2 for the mode structures 



and 1024 2 for the spectrum calculations. Our code is massively parallelized using 
Message Passing Interface (MPI) libraries to facilitate higher resolution. While 
the ISM turbulence code is evolved with time steps resolved self-consistently 
by the coupled fluid motions, the nonlinear interaction time scales associated 
with the plasma 1/k • U«(k) and the neutral 1/k • V n (k) fluids can obviously 
be disparate. Accordingly, turbulent transport of energy in the plasma and the 
neutral ISM fluids takes place distinctively on separate time scales. Small scale 
turbulent fluctuations in ISM evolve under the action of nonlinear interactions 
as well as charge exchange sources. Energy cascades amongst turbulent eddies 
of various scale sizes and between the plasma and the neutral fluids. Because 
of the discrepant nonlinear time scales associated with the plasma and neutral 
fluids, mode structures can be different in the two fluids when they are evolved 
together and in isolation (i.e. decoupled). When the plasma and neutral compo- 
nents are decoupled, the plasma fluid evolves in accordance with the ideal MHD 
where current sheets are typically formed in the magnetic field structures. In the 
presence of charge exchange, the two fluids evolve in a self-consistent manner 
and modify the dynamical properties of both the fluids. One of the most notice- 
able features to emerge from the coupled evolution is that the spectral cascade 
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Figure 1 . Self-consistent simulations of turbulence in the partially ionized 
ISM. Figure shows evolution of the plasma density, velocity magnetic field, 
temperature and pressure. Numerical resolution is 512 2 . 



rates are enhanced substantially in contrast to that of the decoupled case. The 
rapid spectral transfer of energy amongst various Fourier modes in the coupled 
plasma neutral gas leads to smearing off the current sheets in the plasma on a 
shorter time scales. The latter evolves on a relatively long time scales and the 
sheets are formed eventually in the magnetic field. The turbulent equipartition 
is set up in the plasma and the neutral fluid modes. This is shown in Fig 1 
where all the quantities associated with the plasma evolution in our simulations 
are depicted. Seemingly, the neutral fluid, under the action of charge exchange 
sources, tends to enhance the cascades rates by isotropizing the ISM turbulence 
on rapid time scales. One of the plausible reasons of the fast-cascade is that 
the energy transfer rates are enhanced effectively in the coupled ISM turbulence 
system, and that they now require a smaller time to mix the plasma and the 
neutral fluids. In any event, the small-scale sheet-like structures in magnetic 
field compresses (or pinches) the plasma density. Hence the density fluctuations 
develop identically the (thinner than) sheet-like structures that co-exist with 
small-scale turbulent fluctuations in its spectrum as shown in Fig 1. The neu- 
tral fluid, on the other hand, evolves isotropically as stated above by forming 
relatively large-scale structures as shown in Fig 2. 
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Figure 2. Evolution of neutral fluid in the coupled plasma neutral ISM 
turbulence. Shown are the neutral density, velocity, temperature and pressure 
fluctuations. 



3. Energy Cascade Rates in the Coupled ISM Plasma Neutral Gas 

As mentioned above, the plasma neutral fluid coupling in the ISM fluctuations 
introduces charge exchange modes, k ce , that are distinctively different from the 
characteristic turbulent mode k. Typically, k ce /k < 1 in the local ISM. Ac- 
cordingly, the energy cascade time scales are modified as discussed below. The 
nonlinear interaction time-scale in ordinary turbulence is given by 



Tnl 



V(. 



(kv k ) 



where v/~ or V£ is the velocity of turbulent eddies. In the presence of charge 
exchange interaction, the ordinary nonlinear interaction time-scales of fluid tur- 
bulence is modified by a factor k ce /k such that the new nonlinear interaction 
time-scale of ISM turbulence is now 



TNL 



k C e I 

k kvk 



On using the fact that k ce is typically smaller than k, i.e. k ce /k < 1 in the ISM, 
the new nonlinear time is k ce /k times smaller than the old nonlinear time i.e. 
t nl ~ (k ce /k)r n i. This reduced nonlinear interaction time in ISM turbulence is 
likely to enhance turbulent energy cascade rates that are determined typically 
by E^/tnl, where E\. is energy per unit mode. It is beacuse of this reduced 
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interaction time that a rapid spectral transfer of turbulent modes tends to smear 
off the current sheets in the magnetic field fluctuations. 

It is interesting to note that the enhanced cascade rates of ISM turbulent 
modes tend to steepen the inertial range turbulent spectra in both plasma and 
neutral fluids. By extending above phenomenological analysis, one can deduce 
exact (analytic) spectral indices of the inertial range decaying turbulent spectra, 
as follows. The new nonlinear interaction time-scale of ISM turbulence can be 
rearranged as t nl ~ k ce v k /(kv k )(l/kv k ) ~ (r^/r ce ), where r ce ~ {k ce v k )~ l rep- 
resents charge exchange time scale. The energy dissipation rate associated with 
the coupled ISM plasma-neutral system can be determined from e ~ E k /rNL 
relation, which leads to e ~ vf,/(k ce /k 2 v k ) ~ k 2 v\jk ce . According to the Kol- 
mogorov theory, the spectral cascades are local in fe-space and the inertial range 
energy spectrum depends upon the energy dissipation rates and the characteris- 
tic turbulent modes, such that E k ~ e 7 £r . Upon substitution of above quantities 
and equating the power of identical bases, one obtains 

E k ~ e 2 / 3 k- 7 ^ 

plasma spectrum. Similar arguments in the context of neutral fluids, when 
coupled with the plasma fluid in ISM, lead to the energy dissipation rates 
e ~ k 2 v 2 /(k ce /k 2 v k ). This further yields the forward cascade (neutral) energy 
spectrum 

E k ~ e 2 ^k- ll '\ 



4. Conclusion 

In conclusion, one of the most important points to emerge from our studies is 
that the charge exchange modes modify the ISM turbulence cascades dramati- 
cally by expediting nonlinear interaction time-scales, consequently, the energy 
cascade rates are enhanced by isotropizing the ISM turbulence on rapid time 
scales. This tends to modify the characteristics of ISM turbulence which can be 
significantly different from the Kolmogorov phenomenology of fully developed 
turbulence. As shown above, the charge exchange interactions lead to a steeper 
power spectrum. Although turbulence in 2D and 3D possesses distinct spectral 
features characterized essentially by the number of the inviscid quadratic invari- 
ants, our 2D simulations provide a hint that the ISM turbulent spectra in 3D 
case will exhibit a non Kolmogorov-like characteristic. It is to be noted that the 
present model does not consider driving mechanism, hence turbulence is freely 
decaying. Driven turbulence, such as due to large scale forcing, supernova explo- 
tion and CME ejection may force turbulence at larger scales. This can modify 
the cascade dynamics in a manner usually described by dual cascade process. 
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